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The linearized Debye-Hiickel theory for liquid state is shown to provide thermodynamically con- 
sistent virial and energy routes for any potential and for any dimensionality. The importance of this 
result for bounded potentials is discussed. 

Integral equations of liquid theory always involve some approximate closure. 1,2 This is an approximate relation 
between the pair and the direct correlation functions in addition to the exact Orstein-Zernike integral equation. 
Unlike exact theories, these approximate closures introduce well known inconsistencies among different routes to the 
equation of state. Given the knowledge of the pair correlation function, there are clearly many possible routes leading 
to the equation of state, but the most frequently used are the energy, the virial (or pressure), and the compressibility 
routes. 

The degree of inconsistency clearly depends upon the goodness of the approximate closure, so that some closures 
might display weaker differences than others and might even, under some particular circumstances, give no difference at 
all between two particular routes. The most notable example of this (although rarely mentioned in the literature) is the 
virial-energy consistency within the hypernetted-chain (HNC) approximation. 3 Other more recent examples include 
the energy and virial routes in the hard-sphere limit of the square-shoulder potential (for any approximation) , ' and 
again the energy- virial consistency for soft-potentials within the mean-spherical approximation (MSA). 6 

The aim of this Communication is to add one more case to this relatively short list by showing that energy and 
virial routes are completely equivalent within the linearized Dcbye-Hiickel (LDH) approximation for any potential in 
any dimensionality. Our interest in this problem has been triggered by recent investigations on bounded potentials, 7-9 
where this consistency is of particular importance, as we shall discuss. 

Consider an arbitrary potential 0(r) for a homogeneous fluid of TV particles in d dimensions. Newton's third law 
of motion implies that 4>(r) = <f>{— r), but the potential need not be spherically symmetric. The virial equation is 
associated with the compressibility factor Z(p, (3) as 1 

= 1 + ^ J drj,(r)r.V/(r), (1) 

where (3 = 1/ksT is the inverse temperature, p and P are the density and the pressure, respectively, and /(r;/3) = 
e -P<t>(n _ 1 i s the Mayer function. In Eq. (1) we have also introduced the cavity function y(r; p, 0), which is related 
to the pair correlation function g(r; p, (3) by the relation y(r) = e^^^gir). The energy equation is associated with 
the energy U per particle u(p,/3), 
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A standard thermodynamic identity 10 provides the consistency condition between the pressure and energy routes: 

'du\ (dZ_ 



(2) 



d P J \d(3) p (3) 

The explicit notation of the variable being kept fixed and the parameter dependence will be dropped henceforth for 
notational simplicity. It proves convenient to further introduce the quantity w(r; p, (3) = y(r; p, f3) — 1, which is related 
to the potential of mean force ifj(r; p, (3) = —f3~ 1 \ng(r;p,/3) by ~/3[tp(r) — <j>(r)] = ln[l + w(r)]. We anticipate that 
w(r), introduced as a definition at this stage, will assume a particular physical meaning within the LDH approximation 
later on. Then, identity (3) translates into the following condition: 
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Introducing the Fourier transforms w(k; p, (3) and /(k; P) of w(r; p, P) and /(r; /3), respectively, Eq. (4) becomes, after 
standard manipulations, 
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Here we have used w (— k) = w (k) from the symmetry relation <f>(— r) 
convenient form by taking into account the mathematical identity 
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r). Equation (5) can be recast into a more 
df(k)~ 
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Upon integration over k in Eq. (5) the second surface term of the right-hand side of Eq. (6) can be dropped and hence 
we find 
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We remark that no approximations have been carried out so far, and that Eq. (7) is completely equivalent to the 
consistency condition (3). Therefore, any w(k) satisfying Eq. (7) gives thermodynamically consistent results via the 
energy and virial routes. 

We now show that this is in fact the case for the LDH theory which is defined by w(r) = y(r) — 1 where w(k) 
satisfies the scaling relation 11 
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with F (z) = z 2 j '(1 — z). This immediately provides the following expressions 
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Equations (10) and (11) readily yield 
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so that the second integral on the right-hand side of (7) vanishes identically. In addition, Eqs. (9) and (10) show that 
the remaining terms in Eq. (7) are identical. This closes the proof. 

It is useful to put the present result into some perspectives. The LDH theory can be derived from diagrammatic 
methods 1 by summing all simple chain diagrams to all orders in density p. A mathematical device to do this is to 
formally multiply the Mayer function /(r) by a bookkeeping parameter p, and then let p — > 0, so that the leading 
diagrams to be retained at each order are the simple chain diagrams, which then give the dominant contribution to the 
pair correlation function within this approximation. This procedure is physically justified only for bounded potentials 
where |/(r)| can be made arbitrarily small by increasing the temperature, and hence the virial-energy consistency is 
also representative of the exact behavior of the equation of state, unlike the case of unbounded potentials where this 
is not the case and consistency does not automatically ensure exact results. 12 

Representative examples of bounded potentials, recently discussed in the literature, include Gaussian potentials, 13 
penetrable spheres (PS), 7 and penetrable square- well (PSW). 8,9 These potentials are currently of the greatest inter- 
est both from a practical point of view, as they mimic ultrasoft systems such as suitable mixtures of colloids and 
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FIG. 1: Equation of state of the three-dimensional PSW model, Eq. (13), as obtained from the virial route, Eq. (1), (solid line) 
and from the energy route, Eq. (2), (solid circles) for reduced temperature ksT /t a = 8, well width A/a = 0.5, and energy ratio 
£r/ea = 2. Also shown are the results obtained from MC simulations (triangles). 

polymers, 14 and theoretically, as they are compatible with a phase transition even in one-dimensional systems (see 
for instance discussion in Ref. 9). 

We have explicitly numerically checked the virial-energy equation within the PSW model defined by the potential 8,9 



where a is the particle diameter, A is the width of the well, and e r and e a are two positive constants accounting 
for the repulsive and attractive parts of the potential, respectively. Two particles then attract each other through 
a square-well potential of depth — e a < and width A, but can also interpenetrate each other with an energy cost 
e r > 0. Figure 1 depicts both the virial and the energy equation of state for the PSW model at a representative 
state point, as obtained from a numerical solution within the LDH approximation. As expected, we find complete 
numerical consistency, in agreement with the analytical proof. The compressibility equation of state (not shown) lies 
slightly below the energy-virial curve. Figure 1 also includes Monte Carlo (MC) data obtained for the same system 
and state. 15 We observe that at this relatively high temperature the LDH solution provides an accurate equation of 
state, in agreement with the previous discussion on bounded potentials. 

We close this Communication with a few remarks. The analytical proof presented here is patterned after a similar 
proof on the virial-energy consistency within the MSA for a general class of soft potentials which include bounded 
interactions treated here. 6 Given the close relationship between MSA and LDH for soft potentials, 1 the result presented 
here and in Ref. 6 retrospectively can be cast within a unified framework associated with the existence of a scaling 
form in Fourier space akin to Eq. (8). 

Of different nature appears to be the virial-energy consistency within the HNC closure. This is a direct consequence 
of the existence of an explicit expression for the free energy, pressure, and chemical potential as a result of a single 
approximation, thus increasing internal consistency 3,16,17 

On the other hand, the HNC theory can be alternatively viewed as an approximation to the exact diagrammatic 
expansion of the pair correlation function which retains the complete class of particular diagrams (chains, both simple 
and netted, and bundles) and the virial-energy consistency can be also regarded as a direct consequence of this. 2 ' 18 
As the full expansion including all diagrams is of course consistent, an additional further consequence is that the class 
of diagrams not included within the HNC approximation (the so-called elementary diagrams related to the bridge 
function) must also be consistent from the virial-energy point of view. Our result builds upon this argument by adding 
the additional piece of information that the full inclusion of simple chain diagrams only also leads to virial-energy 
consistency. A profound consequence of our result is therefore that the virial-energy consistency is deeply tied to the 
retention of all diagrams within a given class. 19 
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